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We geometrically characterize one-qubit dissipators of a Lindblad type. An efficient parametriza- 
tion in terms of 6 linear parameters opens the way to various optimizations with local dissipation. As 
an example, we study maximal steady-state singlet fraction that can be achieved with an arbitrary 
local dissipation and two qubit Hamiltonian. We show that this singlet fraction has a discontinuity 
as one moves from unital to non-unital dissipators and demonstrate that the largest possible singlet 
fraction is ~ 0.654. This means that for systems with a sufficiently entangled ground state there is 
a fundamental quantum limit to the lowest attainable energy. With local dissipation one is unable 
to cool the system below some limiting non-zero temperature. 


I. INTRODUCTION 

Quantum information theory aims at finding proto¬ 
cols that outperform best classical ones pQ. This can be 
achieved either on a case-by-case basis, studying specific 
problems and finding efficient algorithms, or, one can try 
to characterize the power of a certain generic resource, 
like e.g. entanglement, local interaction etc.. The lat¬ 
ter method can be very powerful, however, if the setting 
considered is too broad, the results can be, from a prac¬ 
tical point of view, less than optimal. In this work we 
study capabilities of an important fundamental resource, 
namely that of local dissipation. 

The most general transformations of quantum states 
are linear completely positive trace-preserving maps 
(CPM) also called quantum channels [2]. A simpler sub¬ 
set consists of transformations that are solutions of local 
Markovian master equations of the Lindblad type [ 211 ]. 
Lindblad equations, whose generic properties have been 
studied in e.g. Eli, are increasingly recognized as a use¬ 
ful resource, be it for universal computation [7], steady- 
state manipulation [8], as an experimental setting, or as 
a tool for nonequilibrium physics. Characterizing the 
set of steady states reachable with given resources is a 
formidable problem with results being usually limited to 
special cases. For pure steady states simple conditions 
determine their stationarity shu, while much less is 
known about mixed steady states, see though Ref. [T2] . 
Frequently an important constraint is locality of inter¬ 
action, either because it is least costly to implement, or 
because it naturally arises in weakly coupled systems. Its 
role on pure steady states has been studied in Ref. m, 
and local translational conservation laws in M- 

In this work we study local one-qubit Lindblad dissipa¬ 
tion. First, we efficiently characterize it, leading to a sim¬ 
ple geometrical picture akin to the celebrated tetrahedron 
geometry of qubit channels m ■ We then demonstrate its 
usefulness by studying the set of states reachable by ar¬ 
bitrary two qubit Hamiltonian and one-qubit dissipation. 
We in particular show that the overlap of the steady state 
with a maximally entangled state, i.e. the singlet frac¬ 
tion, is upper bounded by (3+Vb)/8 « 0.654. This result 
sheds light on the influence of local dissipation on non¬ 


local quantum properties. An important consequence is 
that, provided the system’s ground state is sufficiently 
entangled, local dissipation can not cool the system down 
to low temperatures. There is a “temperature gap” be¬ 
low which one can not go. For the particular setting 
studied this improves on the 3rd law of thermodynam¬ 
ics and various generic zero-temperature unattainability 
results |H H7] or, on ground-state cooling limitations in 
specific situations, e.g. m- We show that not only is the 
zero temperature unattainable, but also that finite (low) 
non-zero temperatures can not be reached with local dis¬ 
sipation. 

II. LOCAL DISSIPATION 

The Lindblad equation is m m 

^ = c(p) = i[ P ,H] + /: Ais (p), (i) 

where £ d i s (p) = 2 LjpL t - pL^Lj - L^Ljp is a su¬ 

peroperator called a dissipator that depends on a set of 
traceless Lindblad operators Lj. The Lindblad equa¬ 
tion, that generates a CPM map via A r = e £r , can 
also be written in a non-diagonal form with a dissipator 

4-fiP) = T.j,k 9 j,kMjpLl-pl]L k -i]i k p), where L, 
is a complete set of orthogonal traceless operators. Ma¬ 
trix gj is called a Gorini-Kossakowski-Sudarshan (GKS) 
matrix [3] and has to satisfy g > 0. 

Expanding Hermitian density operator p on n qubits 
in an orthogonal Hermitian traceless operator basis Fj , 
p = ^ ^j c jFj , so that p is parametrized by real 

coherence vector c, dissipator jCdis(p) = c'j Fj induces 
an affine map c' = Me + t, with real matrix M and 
real vector t. The unitary part i[p, H] of the Lindblad 
equation induces map c' = Nc, with N being real anti¬ 
symmetric. We shall discuss generators £dis instead of 
induced channels A t because of simpler relations and be¬ 
cause they form a convex set while Lindblad channels A r 
do not, e.g. m • 

We are interested in one-qubit dissipators, for which 
M is always symmetric [0], and therefore diagonalizable. 
Therefore, in an appropriate basis £dis can be written in 
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a canonical diagonal form (basis {cr x , cr y , cr z , 1 }), 

{0 0 h\ 

Ad) = 0 -2i±2a 0 t 2 

ciis 0 0 _21 ±32 t 3 ' 

0 0 0 0 / 


( 2 ) 


In its canonical form £// a Lindblad dissipator can be 
parametrized by 6 real parameters, qj and t = (£i, £ 2 , £3). 
The central question is, for which values of these param¬ 
eters ([ 2 ]) is the resulting of a Lindblad form £dis |l]) , 
ie., generates a dynamical semigroup? It is important to 
note - and this is the main advantage of the parametriza- 
tion in Eq. ©> - that £dis is linear in these parameters, 
while on the other hand it is quadratic in Lindblad op¬ 
erators Lj. Having linear parametrization of £di s will 
greatly simplify all optimization problems, including our 
application later on. 

We write the canonical dissipator (I2J) in terms of 
a GKS matrix (written in the basis {cr x , cr y , a z }) 


l ( Q 1 -i*3 ^2 \ 

£ = - hs #2 -hi . ( 3 ) 

\ -h 2 hi #3 ) 

Sufficient and necessary condition for to be of a Lind¬ 
blad form is that g > 0. If the dissipator is unital, t = 0, 
the condition g > 0 translates to three simple conditions 
#1,2,3 > 0. We are now going to show that for non-unital 
case a single additional condition is necessary. 


is not possible to have A = t 2 and Q 7 ^ B because this 
would imply less than 3 positive roots. 

Q — B > 0, together with qj > 0, implies that 
Q > Q - qrt\ ~ # 3^3 > # 1 ^ 1 , in turn leading to q 2 q 3 > £?, 
provided that #1 7 ^ 0. Two similar inequalities hold 
for other index combinations. Summing them together 
means that, if all qj > 0 , Q > B implies A > t 2 . 
For Q = 0 , when one or two qj are zero, inequality 
Q > B must be understood correctly in order to implicate 
A > t 2 : if qi = 0, Q > B directly implies that £ 2,3 = 0 
and then, dividing inequality by gi, we get q 2 q 3 > £?; 
if # 2,3 = 0 all tj = 0. These special cases can be ac¬ 
counted for by instead of Q > B writing the inequality 

t 2 t 2 t 2 

as 1 > — 1 —|- 2 —|-— and understanding that, if any 

— <? 2<?3 <? 1<?2 ° 5 J 

of the denominators is zero, the corresponding tj must 
also b e ze ro. We have proved that if g > 0 then condi¬ 
tions © 5 } hold. In the opposite direction, if qj > 0 and 
Q > B, we have seen that A > t 2 as well as C > 0 and 
all the eigenvalues of g are non-negative. □ 


We see that in addition to simple unital conditions © 
only one additional inequality (J5| is needed to charac¬ 
terize Lindblad dissipators. The situation is similar to 
the one for quantum channels [15 where also a single 
additional condition is required m- 

Important special cases of . - Let us have a look at 
cases when some qj are zero. There are only two possibil¬ 
ities that also exhaust all possible with one Lindblad 
operator (ie., g of rank 1 ): 


Theorem 1 . One-qubit dissipator in the canonical form 
given by Eq. © represents a Lindblad dissipator iff 


Qi > 0, (4) 

y-2 y-2 y-2 

+ ^ + (5) 

#2 #3 # 1#3 # 1#2 


(If any of the denominators in the last equation is zero, it 
must be understood that the corresponding numerator tj 
must also be zero, and the term is left out on the RHS.) 


Proof, g > 0 iff all eigenvalues A j of g are non-negative. 
The characteristic polynomial of Eq. ([3| is p( A) := 
det (g - Al) = 512A 3 - 64CA 2 + 8{A - t 1 ) A - (Q - B), 
where A := qiq 2 + qiq$ + g 2 # 3 , B 1 = qpt\ + q 2 t\ + q 3 tj, 
c := qi + q 2 + # 3 ? Q : = # 1 # 2#3 and t 2 := t • t. Using 
Descartes’ rule of signs we can infer the maximal num¬ 
ber of positive and negative roots by the number of sign 
changes of coefficients in p(X) and p{— A), respectively. 
Because g is Hermitian we known that all roots of p( A) 
are real and we can in fact determine the exact number 
of negative and positive roots. First, for g > 0 all diag¬ 
onal matrix elements qj must be non-negative, qj > 0 , 
which also implies that C > 0 (as tr(g) = (7/8, C = 0 
would imply that = 0). The coefficient in front of A 3 
is positive and in order to have three non-negative roots 
one must also have A > t 2 and Q > B. g has one zero 
eigenvalue if Q = B, if in addition A = t 2 , it has two. It 


(i) Exactly one qj is zero, say q\ = 0. As we have 
seen, this implies t 2 = £3 = 0 and q 2 q 3 > t\. One, 
or if also q 2 q 3 = t\ two, eigenvalues of g are zero. 
An example of such dissipator would be one with a 
single Lindblad operator L = cr x + ia y . 

(ii) Exactly two qj are zero, say q 2 = q 3 = 0. In this 
case all tj must be zero and the dissipator is unital. 
Two eigenvalues of the GKS matrix are zero and, 
up to unitary rotations, we have single Lindblad 
operator L oc cr x . 

If all qj > 0 then one can use inequalities like qiq 2 > t\ 
to show that one can not have A = t 2 (while Q = B, 
ie., g of rank 2, is still possible). The case (ii) is the 

only possibility of with two-times degenerate steady 
state. 


A. Geometry of one-qubit dissipators 


Let us now compare the set of channels obtained from 
Lindblad dissipators to the set of all CPMs. Any one- 
qubit CPM A can be brought to a “diagonal” form [ 2 j, 


A = 


D v 
0 1 


D = diag(Ai, A 2 , A 3 ). (6) 
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FIG. 1. (Color online) The set of A j ([ 6 ]) for one-qubit uni- 
tal Lindbladian channels (colored object). Edges corre¬ 
spond to dissipators composed of one Lindblad operator, 
the surface of two, and the interior from three Lindblad op¬ 
erators. Stripes at the surface, aiding in visualization, are 
isolines of planes perpendicular to vector (1,1,1). Right top: 
crosssection at A 3 = 0.3; full rectangle is the CPM’ tetrahe¬ 
dron boundary. Right bottom: crosssection perpendicular to 
(1,1,1) containing point 0.4(1,1,1); full triangle is the tetra¬ 
hedron’s edge. 


We shall compare the two sets in terms of allowed A j 
(Figs. 0 and [ 2 | and in terms of Vj (Appendix |B|. Con¬ 
ditions which A j and Vj have to satisfy for A to be a 
CPM are well known [T5] . In the unital case (vj = 0) 
A j must lie within a tetrahedron defined by 4 corners 
at A j = ±1 and A 1 A 2 A 3 = 1, whereas for non-unital 
channels an additional inequality has to hold [21]. For 
diagonal Lindblad generator in Eq. ([ 2 ]) it is simple to 

obtain the corresponding Lindblad channel A r = e^dis r 
(without loss of generality we set r = 1) with parameters 
in Eq.0 being given by qi = In h = -WiFxp and 
analogously for other components. For unital Lindblad 
channels the set of all in the space of A j is shown in 
Fig. [T] Because qj > 0 implies A j >0, in the space of A j 
the channel Ai is always in the first octant. By trivially 
concatenating such Lindblad channel by a rotation by 7 r 
around one of the three axes one can change the sign of 
two A j, obtaining symmetrical objects shown in Figs, [l] 
and H The set of unital Lindblad channels is bounded 
by hyperbolic paraboloid surfaces like Ai > A 2 A 3 due to 
qj >0, and is of course smaller than the set of all CPMs 
(tetrahedron). In the 1st octant Lindblad channels fill 
\ of the volume of all CPMs, in the whole tetrahedron 
this fraction is |. For non-unital channels there are three 
additional “shift” parameters Vj. At fixed Vj the set is 
a pinched rounded tetrahedron, the rounding being es¬ 
sentially determined by Eq.([5]), pinching by Eq.Q, with 
logarithms in relations between qj and A j adding addi¬ 
tional complexity. The set of such non-unital Lindblad 
channels, together with a “rounded” tetrahedron set of 
non-unital CPMs (see Ref. m for equations), is shown 


FIG. 2. (Color online) The sets of one-qubit non-unital Lind¬ 
blad channels (solid object) and all CPMs (transparent outer 
shell), both for fixed v =s | (sin 6 cos 0 , sin 0 sin 0 , cos 0) with 
0 — 7t/4 and 0 = 7t/3. The surface corresponds to dissipators 
made from 2, the interior from 3 Lindblad operators. Right: 
same crossections as in Fig. ^ Dotted lines are tetrahedron 
edges (unital CPMs), full curves are crossection with the non- 
unital CPM’s surface (transparent shell), colored region is the 
Lindblad set. 


in Fig. [2] One could also plot the two sets for fixed qj in¬ 
stead 01 for fixed Vj , getting an ellipsoid for the Lindblad 
case, Eq.([5|, for details see Appendix |B| It is interesting 
to note that various tetrahedron-like geometric objects 
often pop up when studying different quantum informa¬ 
tion objects, e.g. [2l H5} l22 ] . a common denominator being 
some form of positivity constraint. 


III. APPLICATION 

The singlet fraction of a 2-qubit steady state poo, 
jC(poo) = is defined as F ■= ('ip\p OQ \ip}^ where | !p) 
(01) — |10))/\/2- States with high singlet fraction can 
be used for quantum processing, for instance teleporta¬ 
tion [23] . We want to find the maximal singlet fraction F 
for a 2-qubit system, maximized over all possible 2-qubit 
H and 1-qubit dissipators £dis (this is different than 
in m where the set is studied for a fixed dissipation). 
The following simple theorem, holding for any number of 
qubits n, will be of help. If = T_1 + CjFj is a 
steady state of Lindblad equation C with dissipator £dis 
and unitary part given by iL, acting as Me + t, then 
Poo := ^rl + k^jCjFj is a steady state of Lindblad 
equation C! having the same H and M as £, while the 
shift in the dissipator C' dis is t' = k t (provided, of course, 
that C dis is also a valid Lindblad dissipator). This can 
be proved by simply verifying that = 0. The 

length of steady-state coherences, i.e., its purity, can be 
increased by increasing t. An important consequence is 
that for fixed M the largest steady-state overlap with 
any given state is reached at the maximal possible shift 
|t|. For 1-qubit £dis the maximal F is therefore reached 
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either for t = 0 (unital dissipator) or for g of rank < 2 . 

If £dis is unital the stationary state coherences c satisfy 
(M + N)c = 0. The steady-state subspace is therefore 
spanned by 1 as well as optional pure states correspond¬ 
ing to nontrivial solutions of (M + N)c = 0. It is well 
known [MU that pure state |y) can be a steady state 
of Lindblad equation iff it is an eigenvector of all Lj as 
well as of H — i JA LjLj. Using Schmidt decomposition 

lx) = E^K)aIA)b we see that L j\ a k) A = S j\ a k) A 

must hold for all k. If the rank of |y) is maximal (for a 
given bipartition such that Lj acts on A) Lj must be an 
identity operator and therefore £dis is zero. In particu¬ 
lar, for a 1 -qubit dissipator pure steady state is always 
separable and therefore F can be at most 1 / 2 . This max¬ 
imum can be reached only if there is a single Lindblad 
operator (two different Pauli matrices do not have a com¬ 
mon eigenvector), e.g. L ~ <j x , and therefore two qj are 
zero. On the other hand, for nonzero shift length t := |t| 
at most one qj can be zero. As a consequence, in the 
limit of small nonzero shift t e the optimal singlet 
fraction can be shown to be 1/4. There is a discontinu¬ 
ous transition in the maximal F from 1/2 to 1/4 as one 
smoothly moves from unital (t = 0 ) to non-unital dissi- 
pators (t 7 ^ 0 ). 

The singlet fraction F is invariant to any local ro¬ 
tation U ® U (U is a 1-qubit unitary) and we can al¬ 
ways rotate dissipator to a basis in which the shift vector 
has only one non-zero component, say t\ ^ 0 ,^ 2,3 = 0. 
Due to symmetry reasons, in the optimal case, we ex¬ 
pect in this basis to have q<± = q% and possibly dif¬ 
ferent q±. Because t also has to be maximal we have 
in addition q^qz = t\ = 1 (because H is arbitrary we 
can set t — 1 without loss of generality). One can ar¬ 
gue (Appendix [a| that for such (acting on the 1st 
qubit) the optimal steady state will be of form p^ = 

fill 2 +Cicr x (J 2 +C 4 I 1 O -2 + C§g\g\ + C\\g\g\ +Ci 3 cr x l 2 . 

One could try to maximize F subject to necessary con¬ 
ditions [12| tr(p^ 0 £^g(p 00 )) = 0, r = 1,2,3, we found it 
though more convenient (Appendix [A| to use an alterna¬ 
tive approach. Operator equation C\poo) = 0 represents 
a set of 15 equations that are linear in parameters dj of 
H = J2j djFj as we U as in c j • They can be written as a 
matrix equation Gd = f , where the matrix G as well as 
the inhomogeneous part f depend on c (but not on d). 
The equation has a solution for d only if f is orthogonal 
to the kernel of G. The kernel of G can be explicitly 
calculated, resulting in three constraints (Appendix |A| . 
The fidelity F = f + c\ + cq + cn, subject to these con¬ 
straints, can be analytically maximized. For fixed ^ 2,3 
the optimum is always at qi = 0, i.e., for single Lindblad 
operator. The absolute maximum is at #2 = Q 3 = 1 when 
F max = (1 + (/?)/4 « 0.6545, where p = (1 + v / 5)/2 is 
the golden mean. Steady state p^ in this optimal case 
is of rank 4 and is entangled. Such optimal p^ can not 
be obtained exactly at ^2 = 1, but only in the limit of 
q2 —y 1, which, however, poses no serious obstacles, see 
Appendix [Aj As we see, local non-unital dissipation can 


create mixed entangled steady-states, even though it can 
not create pure entangled steady-states. Without the lo¬ 
cality constraint considered in this work it is, of course, 
always possible to mitigate entanglement destruction by 
specific dissipation, like for instance when a common 
environment is coupled to non-interacting systems [24] , 
or for high-temperature entanglement in driven [25 or 
steady states [26]. 


IV. CONSEQUENCES 

The fact that maximal singlet fraction F max is bounded 
away from 1 has some important implications. For 
Hamiltonians whose non-degenerate ground-state has an 
overlap with a maximally entangled state larger than 
T max , using only local dissipation one can not “cool” the 
systems down to its ground state. For such H there is a 
fundamental limit to local quantum cooling - a minimal 
temperature below which one can not cool. Note that 
in classical setting, using local cooling by e.g. Langevin 
bath, there are no obstacles to the lowest attainable en¬ 
ergy. In the setting studied this improves on various 
unattainability results for certain specific situations [48] . 
Another interesting consequence is that the steady-state 
F can put a constraint on a possible dissipation, even if 
nothing is known about H or £, namely, if F is larger 
than T max , we immediately know that dissipation can 
not be local. Also, the maximal F is very sensitive to 
unitality. For strictly unital local dissipation it is 1/2, 
whereas for infinitesimally weak violation of unitality it 
drops to 1/4. We expect similar results about lowest 
attainable energy to hold also for more than 2-qubit sys¬ 
tems provided dissipation acts locally on only part of a 
system. The reason is that local dissipation has only a 
limited influence on non-local quantum correlations. As 
a simple example, for an n-qubit Heisenberg chain and 
1-qubit dissipation one can show that a non-degenerate 
steady state is always separable and therefore there is 
again a temperature “gap”. 


V. CONCLUSION 

We have efficiently geometrically characterized single¬ 
qubit Lindblad dissipators, opening the door for vari¬ 
ous optimizations involving local dissipation. To demon¬ 
strate its applicability we have calculated the maxi¬ 
mal singlet fraction achievable by one-qubit dissipation, 
showing that it is less than 1, in turn implying a funda¬ 
mental limit to local quantum cooling. 

I would like to thank B. Zunkovic and H. C. F. Lemos 
for discussions and acknowledge grant P1-0044 from the 
Slovenian Research Agency. 
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Appendix A: Maximal singlet fraction for GKS 
matrix with £ 2=6 = 0 

We want to find out the maximal singlet fraction F 
for with q 2 q z = 1, t\ = 1 and t 2 = £3 = 0, for 
which there are two Lindblad operators, L\ = (y+ 
i cl/y/q 2 )/y/S and L 2 = A necessary condition 

for p to be a steady state of Lindblad equation is that 
C r := tr(p^£^ (poo)) are zero f° r r = 1,2,3 [12]. Let us 
parametrize p as p^ = \t\t 2 - Y^j,k=o c ™> <7 i (T 2i where 
m = j + 4fc + 1 and cr^ 0 ’-” 3 denotes {cr x , cr y , cr z , 1}, re¬ 
spectively. For such parametrization F is 

F = j + ci + cq + cn, (Al) 

i.e., it is given by a sum of coefficients in front of 
<Ji<J 2 , <^ 1^2 5 an( f 0 "i 02 - Conditions 6 r = 0 result in equa¬ 
tions that are of (r + l)-th order in unknown coefficients 
Cm. One could use the method of Langrange multipliers 
to solve this constrained optimization problem. However, 
besides practical solvability issues, there is a fundamental 
difficulty that the domain of allowed c’s is not bounded 
with 61,2,3 = 0 (the reason is that acts only on a 
single qubit). As a consequence, for instance, solving the 
resulting Euler-Lagrange equations for optimization of F 
subject to only C\ = 0, gives a solution for which F = \, 
which, as we shall see, is not the correct maximum. To 
make the domain bounded one could add an additional 
constraint, for instance / := tr(p^) < 1. One difficulty 
with using only C r = 0 though would still remain. 6 r = 0 
is necessary and sufficient if the steady state p has non¬ 
degenerate spectrum, but only necessary otherwise. As 
it turns out, the optimal has in our case a degener¬ 
ate spectrum (one eigenvalue is twice degenerate). We 
therefore use a slightly different approach, where though 
conditions C r = 0 will still be used to first infer that a 
number of coefficients c m are zero in the optimal case. 

Maximizing F subject to C\ =0 as well as I = 
(\ + 4 ^m c m) ^ 1 we have a quadratic maximization 
problem that can be solved exactly. First, one can ob¬ 
serve that in C\ some coefficients come only in perfect 
squares, e.g., 2 (q 3 + qi)(cl + cf 0 + cf 4 + c§) + 2 (q 2 + 
</i)( c 3 + Cj + cf 5 + c 2 L1 ). Therefore, if we have a solu- 
tion with nonzero 02,10,14,3,7,15 if is always better, mean¬ 
ing we will have higher F, to set them to zero and in¬ 
stead increase C 6 ,n- In the maximum we will always 
have 02,3,7,10,14,15 = 0. Then, using Lagrange multipli¬ 
ers one can also show that 05,3,9,12 must as well be zero: 
one has a homogeneous set of linear equations for these 
coefficients with the only solution being 05,3,9,12 = 0, un¬ 
less a Langrange multiplier takes a special value in which 
case though equations for 01,4 do not have a solution. 
Therefore, in the optimal situation only five 01,4,6,11,13 
are nonzero. The fact that only 01,4,6,11,13 appear in the 
steady state with optimal singlet fraction is actually not 
very surprising. F depends on ci, 6 ,n so these coefficients 
will likely appear in the optimal p. Then, dissipator 


couples lil2 to crft 2 as well as licdf to Coef¬ 

ficients 01,4,6,11,13 therefore represent in a way a “mini- 
mar’ set that can satisfy all constraints. Incorporating 
constraints C 2 and 63 into analytical argument is harder 
so we rather show results of numerical optimization. In 
Fig. 3 ] we show the dependence of optimal F on the norm 
ll c A I- 4 Em£{2,3,5, 7 , 8 , 9 , io, 12,14,15} c m- We can see that 
the maximum is reached for 11 c a 11 = 0 , as it was already 
the case for the single constraint C\ = 0 . Adding con¬ 
straints C 2 = 0 and C3 = 0, as well as having q 2 ^ g 3 and 
qi 7^ 0, therefore does not change the conclusion that in 
the optimum only 01,4,6,11,13 are nonzero. Observe also 
from Fig.[ 3 ]that adding condition C3 = 0 to Ci,2 = 0 adds 
very little, for instance, we can analytically show that for 
11 c a 11 = 0 one gets F op t ~ 0.65496 when Ci,2 = 0,1 < 1 , 
while F opt « 0.65451 when 61,2,3 = 0 ,I< 1 . 



FIG. 3. Optimal singlet fraction F for different fixed norm 
11c a 11 2 - Circles show the case with constraints Ci = 0,62 = 
0 ,1 < 1 , and the full curve (almost overlapping with circles) 
shows the case for C\ = 0 , C 2 = 0 , C 3 = 0 ,1 < 1 ; q 2 = q 3 = 1, 
qi 0. 

We are therefore left with 5 nonzero coefficients 
ci,4,6,11,13- Stationary Lindblad equation C(poo) = 0 can 
be written as a matrix equation G(c)d = f(c), for un¬ 
known parameters d of the Hamitonian H = ^ZjdjFj. 
The equation has a solution provided f(c) is orthogonal 
to the kernel of G(c). These give sufficient and neces¬ 
sary conditions on coefficients c in order that the corre¬ 
sponding poo is a steady state. Because we have only 5 
remaining unknown c’s G is sufficiently simple so that 
its kernel, being in general of size 3 , can be analytically 
calculated. First, for fixed q 2 ,3 in optimum one always 
has qi = 0, ie. rank 1 GKS matrix with one Lindblad 
operator. The three kernel conditions are in this case 
<723 Ci - 4 c 4 = 0 , 1 + 2 q 23 ci 3 + D(c% - cf 4 ) = 0 and q 2 c n + 
D(c 4 ce + C11C13) = 0 , with D := 93c 6 /(c 4 cn + c 6 ci 3 ) 
and (£23 : = Q2 T ^3- These are now sufficiently simple 
so that F can be analytically maximized. We also note 
that kernel conditions are stronger than 61,2,3 = 0 . The 
dependence of F on q 2 is shown in Fig. [ 4 ] Not very 
surprisingly, the maximum is achieved for q 2 = q% = 1 . 
Only at that point is 6^ invariant to rotations about 
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q 2 

FIG. 4. The optimal singlet fraction of the steady state for 
1-qubit £^2 with t 2 = t 3 = qi = 0, t\ = 1, q 3 = l/<? 2 - 


the x-axis, similarly as is the singlet state (for q 2 ^ 1 
the symmetry between y and z is lost). The value of 
the maximal singlet fraction is F max = 3+ 8 ^ , reached 
at q 2 = q 3 = 1, while the values of coefficients are 
Cl = c 4 = -Ci 3 = (5 + \/5)/40, c 6 = Cn = l/(4y/5). 
The optimal point of q 2 = q% = 1 for which F = F max 
is in fact degenerate: kernel of G is in this special case 
larger and G(c)d = f(c) has no solutions. Optimality 
can be reached only in the limit q 2 1, which, because 
dependence of F close to q 2 = 1 is quadratic (see Fig. [4]) 
poses no serious obstacle. One possibility is to take 
Cl = F2( 1 + 2< 72/ k )/( 1 + 92) 2 I C4 = 3? 2 (l + 2g 2 /«)/(l + 
6 = FK 1 + 92)/«)Cii = |(1 + g|)/K,d3 = c 4 , 

and H = (k - 2g|)/(4(l - g|)(l + + a\a$\, 

where k ■= \J\ + + 10(/| + Acfc + resulting in a 

singlet fraction F = q 2 /(^q 2 + q 2 + 1 — k) (which is 
not the optimal one for q 2 / 1, but approaches the one 
for q 2 —» 1). Taking q 2 = 1 — e, in the limit e —)> 0, 

the expressions simplify to F = (crfcr^ + ^ 1 ^ 2)5 

and ci = (5 + a/ 5)/40(1 - e 2 /2), c 4 = (5 + \/5)/40, 
c 6 = (1 - e)/(4y / 5), cn = (1 + e)/(4>/5), c 13 = -c 4 , 

all written to the lowest order in e. 


Appendix B: Comparing Lindblad and quantum 
channels 


v, such that the resulting map is completely positive or 
can be generated by Lindblad evolution, Eq. (5). Maxi¬ 
mal allowed length of v in a given direction is therefore 
equal (upto a scale) to the distance between the shown 
surface and its center, i.e., A, around which the “ball” is 




FIG. 5. (Color online) The sets of allowed channel shifts v 
[Eq. (6)] (scaled by factor 0.3) for 7 different diagonal val¬ 
ues A, namely for (0,0, 0), (0, 0,0.6), (0,0,0.9), (0.5,0.5,0.5), 
(0.8,0.8, 0.8), (0.5, 0,0) and (-0.6,-0.6,0.4). Top: Lindblad 
channels. Bottom: CPMs. What is shown is a polar plot: 
maximal v in a given direction, centered around fixed A. In 
the Lindblad case the surface of “balls” corresponds to dissi- 
pators with 2 Lindblad operators, and the interior to 3. On 
the diagonal A oc (1,1,1) the set of allowed v is the same for 
CPMs and Lindblad channels (perfect spheres), everywhere 
else the Lindblad set is smaller. 


In the main text we have compared the set of single¬ 
qubit Lindblad channels to the set of general qubit chan¬ 
nels for the unital case in Fig. 1, and for a fixed shift vec¬ 
tor v, Eq. (6) in main text, in Fig. 2. It is instructive to 
compare the two sets in a non-unital case also for chan¬ 
nel’s fixed generalized singular values A = (Ai,A 2 ,As), 
Eq. (6). 

These sets are shown in Fig. [5] What is shown is a po- 
lar plot of the maximal allowed length of the shift vector 


plotted. For Lindblad channels, in the space of genera¬ 
tor shifts t, this set would be an ellipsoid, Eq. (5). In 
the space of channel shifts v it is a smooth ellipsoid-like 
shape seen in the top of Fig. [5j For general channels 
(CPMs) the set, which is larger or equal to the Lindblad 
one, can be seen in the bottom of Fig. [5] The boundary 
of this CPM set is determined by Eq.(30) from Ref. [ 24) 
and can have non-smooth edges. 


[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information , (Cambridge, 2000). 


[2] I. Bengtsson and K. Zyczkowski Geometry of Quan- 





















7 


turn States: An Introduction to Quantum Entanglement , 
(Cambridge, 2006). 

[3] V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, 
Completely positive dynamical semigroups of N-level sys¬ 
tems , J. Math. Phys. 17, 821 (1976). 

[4] G. Lindblad, On the generators of quantum dynamical 
semigroups , Commun. Math. Phys. 48, 119 (1976). 

[5] B. Baumgartner, H. Narnhofer, and W. Thirring, Anal¬ 
ysis of quantum semigroups with GKS-Lindblad genera¬ 
tors: I. Simple generators , J. Phys. A 41, 065201 (2008); 
B. Baumgartner and H. Narnhofer, Analysis of quantum 
semigroups with GKS-Lindblad generators: II. General , 
J. Phys. A 41, 395303 (2008). 

[6] F. Ticozzi and L. Viola, Analysis and synthesis of attrac¬ 
tive quantum Markovian dynamics , Automatica 45, 2002 
(2009). 

[7] F. Verstraete, M. M. Wolf, and J. I. Cirac, Quantum com¬ 
putation and quantum-state engineering driven by dissi¬ 
pation , Nature Phys. 5, 633 (2009). 

[8] P. Zanardi and L. C. Venuti, Coherent quantum dynamics 
in steady-state manifolds of strongly dissipative systems , 
Phys. Rev. Lett. 113, 240406 (2014). 

[9] P. Zanardi, Dissipation and decoherence in a quantum 
register , Phys. Rev. A 57, 3276 (1998). 

[10] N. Yamamoto, Parametrization of the feedback Hamil¬ 
tonian realizing a pure steady state , Phys. Rev. A 72, 
024104 (2005). 

[11] B. Kraus, H. P. Biichler, S. Diehl, A. Kantian, A. Micheli, 
and P. Zoller, Preparation of entangled states by quantum 
Markov processes , Phys. Rev. A 78, 042307 (2008). 

[12] S. Sauer, C. Gneiting, and A. Buchleitner, Optimal co¬ 
herent control to counteract dissipation , Phys. Rev. Lett. 
Ill, 030405 (2013); S. Sauer, C. Gneiting, and A. Buch¬ 
leitner, Stabilizing entanglement in the presence of local 
decay processes , Phys. Rev. A 89, 022327 (2014). 

[13] F. Ticozzi and L. Viola, Stabilizing entangled states 
with quasi-local quantum dynamical semigroups , 
Phil. Trans. R. Soc. A 370, 5259 (2012); F. Ticozzi 
and L. Viola, Steady-state entanglement by engineered 
quasy-local Markovian dissipation , Quantum Information 
and Computation 14, 0265 (2014). 

[14] M. Znidaric, G. Benenti, and G. Casati, Translationally 
invariant conservation laws of local Lindblad equations , 
J. Math. Phys. 55, 021903 (2014). 

[15] M. B. Ruskai, S. Szarek, and E. Werner, An analysis of 
completely-positive trace-preserving maps on 2x2 matri¬ 
ces , Lin. Alg. Appl. 347, 159 (2002). 

[16] L.-A. Wu, D. Segal, and P. Brumer, No-go theorem for 
ground state cooling given initial system-thermal bath 
factorization , Sci. Rep. 3, 1824 (2013); F. Ticozzi and 
L. Viola, Quantum resources for purification and cool¬ 


ing: Fundamental limits and opportunities , Sci. Rep. 4, 
5192 (2014); 

[17] L. Masanes and J. Oppenheim, A derivation (and 
quantification) of the third law of thermodynamics , 

arXiv:1412.3828 (2014). 

[18] T. Feldmann and R. Kosloff, Minimal temperature of 
quantum refrigerators , EPL 89, 20004 (2010); X. Wang, 
S. Vinjanampathy, F. W. Strauch, and K. Jacobs, Ab¬ 
solute dynamical limit to cooling weakly coupled quantum 
systems , Phys. Rev. Lett. 110, 157207 (2013); G. Benenti 
and G. Strini, Dynamical Casimir effect and minimal 
temperature in quantum thermodynamics , Phys. Rev. A 
91, 020502(R) (2015). 

[19] M. M. Wolf, J. Eisert, T. S. Cubitt, and J. I. Cirac, 

Assessing non-Markovian quantum dynamics , 

Phys. Rev. Lett. 101, 150402 (2008). 

[20] K. Lendi, Evolution matrix in a coherence vector formu¬ 
lation for quantum Markovian master equations for N- 
level systems , J. Phys. A 20, 15 (1987). 

[21] D. Braun, O. Giraud, I. Nechita, C. Pellegrini, and 
M. Znidaric, A universal set of qubit quantum channels , 
J. Phys. A 47, 135302 (2014). 

[22] I. Bengtsson, S. Weis and K. Zyczkowski, Geometry of the 

set of mixed quantum states: An apophatic approach , in 
Geometric Methods in Physics. XXX Workshop, pp 175- 
197 (Springer, Basel, 2013); P. D. Johnson and L. Vi¬ 
ola, On state versus channel quantum extension prob¬ 
lems: Exact results for U symmetry , J. Phys. A 

48, 035307 (2015). 

[23] M. Horodecki, P. Horodecki, and R. Horodecki, General 
teleportation channel, singlet fraction, and quasidistilla¬ 
tion , Phys. Rev. A 60, 1888 (1999). 

[24] D. Braun, Creation of entanglement by interaction with 
a common heat bath , Phys. Rev. Lett. 89, 277901 (2002). 

[25] F. Galve, L. A. Pachon, and D. Zueco, Bringing Entan¬ 
glement to the high temperature limit , Phys. Rev. Lett. 
105, 180501 (2010). 

[26] L. Quiroga, F. J. Rodriguez, M. E. Ramirez, and R. Paris, 
Nonequilibrium thermal entanglement , Phys. Rev. A 75, 
032308 (2007); L.-A. Wu and D. Segal, Quantum ef¬ 
fects in thermal conduction: Nonequilibrium quantum 
discord and entanglement , Phys. Rev. A 84, 012319 
(2011); M. Znidaric, Entanglement in stationary nonequi¬ 
librium states at high energies , Phys. Rev. A 85, 012324 
(2012); B. Bellomo, and M. Antezza, Creation and pro¬ 
tection of entanglement in systems out of thermal equi¬ 
librium , New J. Phys. 15, 113052 (2013); S. Camalet, 
Steady Schrodinger cat state of a driven Ising chain , 
Eur. Phys. J. B 86, 176 (2013). 


